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ABSTRACT 

V 

Inteznal  stresses  inside  a prolate  s^dieroidal  inclusion  parallel  to 
a uniaxial  ^iplied  stress  are  obtained  by  using  Eshelby  theory.  Effects 
of  the  elastic  nodulus  and  the  cispect  ratio  of  the  inclusion  are  evaluated. 
Ihe  present  results  are  ocii|)ared  with  others*  in  particular*  with  Argon's 
solution  for  a slender  rod. 
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1.  INTRODUCnCN 

When  exteznal  stress  is  apjplied  to  a oaoposite  material,  internal 
stresses  witMn  the  material  are  strongly  affected  by  tiie  shape,  arranganent, 
and  elastic  stifftiess  of  the  reinforcing  element.  For  (^imum  design  and 
effective  utilization  of  ooa|X3Bite  materials,  Ihorou^  knowledge  of  the 
stress  oonoentraticn,  eq>ecially  at  tiie  interface,  is  essential.  It  is 
pertinent  to  note  that  failure  often  ini-^tes  at  such  stress  oonoentraticn. 

In  fiber  reinforced  oon|X3site  materials,  the  stress  concentration  at  the 
ends  of  disocntinuous  fibers  is  of  prijne  ioqportanoe.  This  problem  Ms  been 
oonsidered  by  several  investigators  Cl*'|].  Most  recently.  Argon  obtained 
e^iproximate  solutix^  for  stresses  around  slender  elastic  rods  with  moduli 
different  from  that  of  the  infinite  elastic  matrix  [1].  Earlier,  Edwards  [2] 
evaluated  the  stress  concentration  due  to  a prolate  ^heroidal  inclusion 
which  lies  normal  to  the  applied  stress  axis.  TanaJca  et  al.  [3]  oonsidered 
a prolate  spheroidal  inclusion  parallel  to  the  stress  a:<is  in  the  limit  of  an 
infinitely  long  inclusion.  While  a prolate  spheroid  is  a good  approximation 
for  a fiber,  no  detailed  study  of  the  stress  oonoentration  has  been  made.  In 
this  stuc^,  we  report  on  the  stress  concentration  produced  by  a prolate 
spheroidal  inclusion,  which  lies  in  an  infinite  elastic  matrix  parallel  to 
the  applied  stress  axis.  Effects  of  the  shape  and  the  elastic  noduli  of  the 
inclusion  on  the  stress  concentration  are  considered.  Eshelby  theory  [4] 
of  transformation  induced  stresses  foams  the  basis  of  the  present  stucty. 

This  ^3proach  is  ooomon  to  our  previous  studies  on  the  stress  concentration 
in  arKl  around  an  oblate  spheroidal  indusion  embedded  in  the  elastic-plastic 
matrix  [5,6]. 

2.  mCLUBION  /HD  /tfUED  STBESS 

In  the  CartMlaA  aoordinebea  eystm,  X2  and  X^,  a prolate  spheroidal 
inclusion  is  loostsd  within 

0^  ♦ J^)/a*  ♦ J^/c*  < 1 . 


(1> 


for  inclusion,  respectively.  XandX*  are  the  Lams  constants,  and  y and  y* 
are  the  shear  npduli  of  the  matrix  and  uiclusion,  respectively.  In  ocder  to 
si]i|)lify  the  presentation  of  rumerioal  z^esults,  Poisson's  ratios  of  the 
matrix  and  inclusion,v  and  v*  are  assumed  to  be  equal  to  1/3,  unless  othezvise 
specified.  In  the  following  calculations , we  employ  the  ratio  of  shear  modulus 
m s y*/y  under  uniaxial  applied  stress  along  the  axis.  The  elastic 
strain  e^.  produced  by  oan  be  written  as. 


where  E is  the  Young's  modulus  of  the  matrix 


3.  THEORY 


Eshelby  [4]  obtained  general  soluticn  for  the  ihhonogenity  problem  of 
ellipsoidal  inclusion  using  the  ao-oalled  equivalent  inclusion  metliDd 
internal  stresa  inside  the  inclusion,  which  is  uniform,  is  given  by 


'ij  * ^ijJdl^*£l  * * ^ijldl^®kA  * 


vAiere  the  constraint  stredn  ^ is  given  in  terms  of  Edielly  tensor, 
c T T • 

as,  ^ *im  “ eigen  strain  of  the  equivalent  inclusion. 

By  solving  (3)  with  respect  to  e^ , vrs  can  obtain  o?, . Because  of  the  symnetxy  | 

m I 


of  Eqs,  (1)  and  (2),  tates  the  fom  of 


w 


vAiere  X and  Z sre  functic^  of  1i»e  aspect  ratio,  k,  and  the  ratio  of  shear 
modulus,  m,  as  well  as  d^/E. 


4.  RESULTS 

From  Eq.  (3),  X and  Z are  found  to  be 


Y - (1  -»  m)  [(1  - a)  r(k)  - IJ 

X -j—— 


3t(l  - jn)^f(k)  - (1  - in)g(k)  + 1] 


Z s -(1  - m)  [(1  - m)  s(k)  - 3]  . / 

3[(1  - m)^f(k)  - (1  - m)g(k)  + 1]  ^ 


(5) 


where  m is  y*/y,  and 

* ^3333  * ^®1133 
8(k)  = + 3Sjj22  ^3311 

f(k)  s ♦ ®1122^^3333  “ ^1133^3311 

g(k)  s S^^  + S^22  * ^3333  * 

The  values  of  r,  s,  f and  g have  been  evaluated  for  k ^ 1 by  using  appropxdate 
^Ijkft  resultant  internal  Ctrssses  inside  the  inclusion  can  be  written 


<^  * <^2  ■ 2y*  . * <3X  ♦ Z) 


. 3 . 


te 


it 


<^3  = “I'*  • rh:  • « * z)  • <6) 

The  nunerioal  results  are  shoNn  in  Figs.  1 and  2.  The  stresses  are 
nomialized  by  and  are  plotted  against  k for  selected  values  of  m. 

Results  indicate  that  both  and  rapidly  approach  the  asvnptotic 
values  vdth  increasing  k and  significant  deviations  are  only  found  when  k 
is  less  than  10  to  30.  For  an  infinitely  long  inclusion  always  vanishes, 
and  becomes  ntA.  When  no  inhooogenaity  effect  exist  (m  = 1),  we  obtain 
, and  ^ 0 as  expected.  The  shear  modulus  of  the  inclusion  has  a 
larger  effect  on  for  a more  elongated  inclusion. 

5.  00rf>AK[S0N  Wim  OTHER  STUDIES 

5-1.  The  Results  for  an  Infinitely  long  Inclusion  (k  = ») 

Tanaka,  Mori  and  Nakanura  made  a pioneering' contribution  to  inclusicxi 
problems  by  using  Eshelby  Iheory.  In  one  of  their  papers  [3],  the  internal 
stress  due  to  Ihe  inhomogeneity  effect  of  an  infinitely  long  fiber  parallel 
to  the  stress  axis  was  evaluated  as, 

= {1  + (E*  - E)  [(1  ♦ V)  (1  - 2v)E  * Cl  - v)E*]  j ^ 

(1  - v^)  C(1  - 2v)E  + E*]E 

Their  result  disagrees  with  the  present  calculation  except  then  specific 
values  of  m are  taken.  Relevant  non-zero  conponents  of  Eshelby  tensor  for 
a prolate  spheroidal  inclusion  parallel  to  for  k s • are 

C a c - 5 - 4v 
^1111  **2222  ‘ 8(1  - V) 

^1122  ^2211  8(1  -'v) 


froB  which  we  have  e. 


5-2.  Rad-Shaped  Fiber 

Recently,  Argcn  [23  obtained  an  f^^praxiioate  solution  for  Ihe  stress 
oonoentrsticn  due  to  a rod-shaped  inclusion  in  an  infinite  elastic  mediun 
when  uniaxial  external  stress  is  applied  along  the  rod-axis.  He  used 
Esheltv's  concept  of  stress-free  transforaation  in  order  to  evaluate  traction 
at  interfaces.  Then,  he  enf>loyed  a solution  of  elasticity  for  two  point 
forces  in  an  infinite  medium  given  by  Timo^vecto  and  Goodier  [8].  The 
interfeKse  traction  was  represented  by  the  two  point  forces  concentrated  at  the 
end  of  the  rod-shaped  inclusion.  A rod  having  a large  length- to-diameter  ratio 
is  siniLar  to  a prolate  spheroidal  inclusion  tdthalarge  aspect  ratio.  There- 
fore, it  is  interesting  to  exasdne  differences  between  results  obtained  by 
Argon  and  by  us. 

When  a rod-shaped  inclusion  with  a diameter  2 r^  and  length  21,  is 
parallel  to  the  axis  of  extemal  stress  cA,  Argon's  solution  for  at  the 
end  of  the  rod  is 


«^)ere  v is  assuniBd  to  be  1/3,  and  a constant  a is  found  to  be  1.412.  Here, 
the  apseot  ratio  k*  is  defined  as  k*  « l/r  . In  Fig.  3,  values  of 
obtained  fxom  Eq.  (6)  and  Eq.  (9)  are  plotted  asdnst  k and  k*  for  m a 6.  (This 
value  of  m oorresponds  to  a steel  reinforced  concrete. ) As  can  be  seen  in 
Fig.  3,  a good  agreement  was  obtained.  When  k (or  k’)  goes  to  infinity,  both 
results  agree  exactly.  Differenoes  are  less  tl^  3%  for  k (or  k* ) > 100  and 
becane  greater  with  decreasing  k (or  k’ ) . However,  it  is  still  less  than 

t 

20%  at  k = 1.  As  pointed  out  by  Argon,  his  method  based  on  two  point  forces 
in  an  infinite  mediun  is  applicable  to  'ttie  Ccise  of  very  large  aspect  ratios 
of  tJie  rod-shaped  incliision.  A part  of  liie  difference  should  eLLso  be  attri- 
buted to  different  geometries. 
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FIGURE  CAPTIONS 


Fig.  1.  Internal  stress  inside  the  inclusion  due  to  inhonogeneity  effect, 
against  k. 

Fig.  2.  Internal  stress  inside  live  inclusion  due  to  inhonogeneity  effect, 
aJa//,  against  k. 

Fig.  3.  Internal  stress  inside  the  prolate  and  rod-shaped  inclusions 

at  m = 6. 
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